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ABSTRACT 

Context. Supermassive black holes (SMBHs) are fundamental keys to understand the formation and evolution of their 
host galaxies. However, the formation and growth of SMBHs are not yet well understood. One of the proposed forma- 
tion scenarios is the growth of SMBHs from seed intermediate-mass black holes (IMBHs, 10^ to 10^ Mq) formed in star 
clusters. In this context, and also with respect to the low mass end of the AI, — a relation for galaxies, globular clusters 
are in a mass range that make them ideal systems to look for IMBHs. Among Galactic star clusters, the massive cluster 
Lj Centauri is a special target due to its central high velocity dispersion and also its multiple stellar populations. 
Aims. We study the central structure and dynamics of the star cluster uj Centauri to examine whether an IMBH is 
necessary to explain the observed velocity dispersion and surface brightness profiles. 

Methods. We perform direct N-body simulations on GPU and GRAPE special purpose computers to follow the dy- 
namical evolution of lj Centauri. The simulations are compared to the most recent data-sets in order to explain the 
present-day conditions of the cluster and to constrain the initial conditions leading to the observed profiles. 
Results. We find that starting from isotropic spherical multi-mass King models and within our canonical assumptions, 
a model with a central IMBH mass of 2% of the cluster stellar mass, i.e. a 5. x lO^M© IMBH, provides a satisfactory fit 
to both the observed shallow cusp in surface brightness and the continuous rise towards the center of the radial velocity 
dispersion profile. In our isotropic spherical models, the predicted proper motion dispersion for the best-fit model is 
the same as the radial velocity dispersion one. 

Conclusions. We conclude that with the presence of a central IMBH in our models, we reproduce consistently the rise in 
the radial velocity dispersion. Furthermore, we always end up with a shallow cusp in the projected surface brightness 
of our model clusters containing an IMBH. In addition, we find that the M/L ratio seems to be constant in the central 
region, and starts to rise slightly from the core radius outwards for all models independent of the presence of a black 
hole. Considering our initial parameter space, it is not possible to explain the observations without a central IMBH 
for uj Centauri. To further strengthen the presence of an IMBH as a unique explanation of the observed light and 
kinematics more detailed analysis such as investigating the contribution of primordial binaries and different anisotropy 
profiles should be studied. 

Key words. Black Holes - globular clusters: individual (Omega Centauri) - stellar dynamics - methods: N-body simu- 
lations 



1. Introduction star clusters and IMBHs in the bulge. They find that 10% 

^ ' of the clusters born within 100 pc of the Galactic center 

d I There is no doubt about the existence of supermassive black undergo core collapse during their inward migration and 

holes (SMBHs) at the center of most galaxies. However, form IMBHs via runaway stellar merging. The IMBHs con- 

the formation and growth of SMBHs is poorly understood, tinue their inward drift towards the Galactic cent er after 

One of the proposed scenarios is the growth of SMBHs from the di ssolution of the host clusters. Portegies Zwar t et all 

seed intermediate-mass black holes (IMBHs , 10^ t o IO^Mq) (|2OO60 predict that a region within 10 pc of the Galactic 

(jEbisuzaki et all 1200X1: iTanaka fc HaimanI I2009D . IMBHs center might be populated by 50 IMBHs of about 1000 Mq 

in star clusters might form through the runaway merg- mass. They also predict that there is a steady population 

ing of massive s tars ("Porte gies Zwart fc McMillanI l2002t of several IMBHs within a few miUiparsecs of the Galactic 

iPortegies Zwart et al. 2004). center. This population merges with a rate of about one 

IMBH formation in star clusters could help to ex- Per 10 Myr with the central SMBH, which is sufficient 
plain the supermassive black hole formation and gr owth in build the accumulated majority of the SMBH mass. In 

the ce nter of galaxies. For instance, Portegies Zwart et al.l t^^e same context nuclear star clusters co-exist with mas- 

(|2006l) simulate the inner 100 pc of the Milky Way to sive black holes (jSeth et al.| [20081, 12010). The star cluster 
study the formation and evolution of the population of Centauri (NGC 5139) in our galaxy might be a bridge 
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between smaller systems such as classical globular clusters 
and larger systems like nuclear star clusters. 

If IMBHs form in large numbers in star clusters, then 
one might expect that some star clusters in the Milky Way 
or other nearby galaxies contain central black holes. For 
example, iLiitzgendorf et all (|201lh report a kinematic ev- 
idence for the exsitence of an IMBH in NGC 6388. The 
structural paramete rs of globular clust e rs harb or ing IM BHs 
are also studied in iBaumeardt et al.l (|2004alibl . 120051 ) and 
IVesperini fc Trentil (|2010t) . Care should be taken when in- 
terpreting cluster mor phological par ameters as IMBH in- 
dicators. For instance, 'Hurlev (2007) shows that the large 
cores observed in some Galactic star clusters can be caused 
by heavy stellar mass black hole binaries without the need 
to invoke an IMBH. Baumgardt et al. (2005) show that 
core-collapsed globular clusters with steep surface bright- 
ness profiles are not good candidates for harboring cen- 
tral black holes. They find that a cluster hosting an IMBH 
appears to have a relatively large core with a projected 
surface brightness only slightly rising toward the center. 
It should be noted that Vesperini & Trenti (2010) ar- 
gue that shallow cusps in the central surface brightness 
profile may not be a uniq ue IMBH indicator (see also 
iNovola fc Baumgardt! (|201lh for a different interpretation). 
Baumgardt et al. (2005) show that the velocity dispersion 
of the visible stars in a globular cluster with a central black 
hole remains nearly constant well inside the apparent core 
radius. Further, they report that in a cluster containing an 
IMBH, the influence of the black hole becomes significant 
only at a fraction | • of the half-mass radius (where 

Mbh and Mc are the mass of the IMBH and the clus- 
ter), i.e. deep within the core, which will affect only a small 
number of stars. 

The star cluster a; Centa uri (NGC 5139), with a n esti- 
mated mass of 2.5 x lO^M© (jvan de Ven et al.l 12006 here- 
after vdV06) and a tidal radius of about 70 pc (Harris 
1996), is the most massive and one of the most spatially ex- 
tended Galactic star clusters. It has one of the highest cen- 
tral velocity dispersi ons among the Milky Way star clusters 
with about 22 km/s (jMevlan et al.lll995l:lNovola et al.l20Tol 
hereafter NIO). Furthermore, vdV06 measure a rotation of 8 
km/s at a radius of about 11 pc from the center using radial 
velocities. In addition, ui Centauri is one of the first Galactic 
globular clusters that have multiple stell ar populations 
among both red giant branch stars ([Freem an & Rodgers 
1975 ) and main sequence stars ()Andersonll2002tlBedin et al 
2004f) . The nature of this cluster is therefore a matter of de- 



bate, it could either be a gia nt globular cluster or the core 
of a stripped dwarf galaxy (jFreemanl 119931 : iMevlanI [2002t 
iBekki fc Norrisl[2b06l) . The above spectacular properties, in 
addition to a shallo w cusp in the surface brightness profile 
()Novola et al.ll2008l hereafter NGB08) and a central sharp 
rise in the radial velocity dispersion (NIO), make uj Centauri 
an interesting candidate for harboring a black hole. 

w Centauri's dynamics is among the best studied of any 
Galactic star cluster. vdV06 determine its dynamical dis- 
tance, inclination, mass-to-light ratio, and the intrinsic or- 
bital structure by fitting axisymmetric dynamical models 
to the ground-based proper motions of van Leeuwen et al. 
(2000) and line-of-sight velocities from independent data- 
sets. They find that w Centauri shows no significant radial 
dependence of M/L, consistent with its relatively long re- 
laxation time. Their best-fit dynamical model has a stellar 
V-band M/L of 2.5 ± 0.1 (solar units) and an inclination 



i — 50° ± 4°, which corresponds to an average intrinsic ax- 
ial ratio of 0.78 ± 0.03. These models do not include an; 
kinematical data in the central 10". Gicrsz & Hcggic (200 
use Monte Carlo simulations to model uj Centauri with sim- 
ple spherical models (neglecting rotation and binary stars) . 
They fit the surface brightness and radial velocity disper- 
sion relatively well, though again, neither the data nor the 
model have sufficient resolution in the central 1 pc (~ 43" 
at 4.8 kpc). 

There are several well established correlations be- 
tween the central black hole mass of galaxies and 
other pa rameters of host galaxies such as velocity dis- 
persion (iFerrarese fc MerrittI [2OOOI : iGebhardt etHI I2OOOI: 
iGiiltekin et al.l I2009D. bulge mass and bu lge luminosity 
(|Magorrian et all 119981: iHaring fc Rixll2004{ ). If we extrap- 
olate the Magorrian et al. (1998) relation to the globu- 
lar cluster mass regime, it predicts an IMBH of about 
1.5 X 1O^M0 for UJ Centauri assuming a total cluster mass 
of 2.5 X IO^Mq. NGB08 find a (4 ± 1) x IO^Mq IMBH 
applying isotropic Jeans models and a (3 ± 1) x 10'* M0 
IMBH using axisymmetric orbit-based models. More re- 
cently, NIO provide new central kinematics of w Centauri 
and suggest a (5.2±0.5) x lO'^M© IMBH assuming spherical 
isotropic Jeans models with res pect to a newly determined 
kinematic center. In contrast, Ivan der Marel fc AndersonI 
(2010. , hereafter vdMAlO), using HST proper motions find 
a (8.7 ± 2.9) x 1O^M0 IMBH assuming cusp models and an 
upper limit of 7.4 x IO^Mq at la confidence assuming core 
models (flat central density), while isotropic models imply 
an IMBH mass of (1.8 ± 0.3) x IO'^Mq. One of the main 
reasons for the discrepancy is the different cluster centers 
these two groups used. In addition to the center determi- 
nation, underestimating the rotation, particularly in the 
central parsec, could have an important effect on velocity 
dispersion measurements. 

In this paper, we compare the most up-to-date observed 
surface brightness and kinematic profiles of uj Centauri with 
direct N-body simulations in the same way as observers 
do. This means the same luminosity weights and magni- 
tude cut-offs as in the observations are applied to compute 
the velocity dispersion and surface density profiles from the 
models. Similar studies have been performed earlier for the 
globular cluster M 15 in the Milky Way and Gl in M31 by 
Baumg ardt et all ([2003a,b). Direct N-body simulations of 
Ml 5 explain the observations with a concentration of dark 
remnants, such as massive white dwarfs and neutron stars 
in the central regions through mass segregation. Therefore, 
the presence of an IMBH was not necessary in order to ex- 
plain the observatio n s. Th e same conclusion is made by 
Ivan den Bosch etHI (l2006l) for M15 using Schwa rzschild 
model. In the case of Gl. iBaumgardt "eTaLl (l2003bD repro- 
duce the observations by a s sumin g a merger history for Gl. 
However, IGebhardt et al.l (|2005|) provide additional sup- 
port for the presence of a 2 x 10^ Mq IMBH. The black 
hole scenario for Gl is also supporte d by detections of ra- 
dio and x-ray sou rces in the cluster (|Ulvestad et al.l [20071 : 
iPoolev fc Rappap ort 2006). 

Here, we examine different IMBH masses (including the 
no black hole case) in our N-body models with the aim of 
reproducing the observations for uj Centauri. Only iV-body 
models allow realistic inclusion of relaxation and stellar evo- 
lution effects and changes in M/L with radius due to mass 
segregation. The possible disadvantage of N-body models 
is that one is restricted to a small number of models since 
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they are time consuming to construct. Therefore, we re- 
strict ourselves to isotropic models and only run a three 
dimensional grid in concentration (c), projected half- mass 
radius {rtp) and IMBH mass space for the scope of this pa- 
per. Although axisymmetry and anisotropy are important 
to include in any modeling, according to vdV06 w Centauri 
is close to isotropic and spherical within the central few 
core radii. 

In Section 2, we describe the data used in this work to 
compare with the N-body results. The general recipe for 
our N-body models is discussed in Section 3. In Section 4 
we explain our model results, in particular we discuss the 
profiles of the no-IMBH and IMBH models and compare 
them with the observations. We draw our conclusions and 
discuss possible future work in Section 5. 

2. Observational Data 

2.1. The Center of uj Centauri 

The determination of u Centauri's center is crucial in order 
to understand and model its dynamics. However, the exact 
location of the center has been controversial due to the 
large flat core of the cluster (core radius ^ 100") and the 
different methods used to estimate its location. 

The center of w Centauri has been determined by sev- 
eral authors. Recent determinations are done by NGB08 
and [Anderson fc van der Marell (|2010l hereafter AvdMlO). 
NGB08 determine the center of w Centauri using star 
counts by excluding the faintest stars due to incomplete- 
ness. This method can be biased towards bright stars. In an 
independent study, AvdMlO determine the center with dif- 
ferent methods including star counts and proper motions. 
In this case, the authors use star lists corrected for the 
presence of bright stars assuming a symmetry axis. This 
measurement can be biased due to the quality of the correc- 
tion and the location of symmetry axis. Their result differs 
from the NGB08 position by ~ 12". AvdMlO also deter- 
mine the center of u Centauri using HST proper motion 
data, which they report to be in agreement with their star 
count method within the uncertainties. AvdMlO might have 
underestimated the rotation contribution in their local fil- 
ter window in the proper motion measurements. They try 
to estimate global rotation but were limited in the amount 
they could detect. Their evaluated center based on proper 
motion could possibly be offset from the true center due to 
this effect. 

Due to the above discrepancy, NIO argue that using the 
kinematic center rather than the density center is the better 
choice as starting point for models. For our N-body models 
we use the kinematical center derived in NIO. 

2.2. Surface Brightness Data 

iMevlanI (|l98iD and iTraeer et all (|1995D compile surface 
brightness data for w Centauri from different sources in 
the li terature: aperture photometry of the central regions 
from iGascoigne fc Buril (|l956h ai id iDa Costal (Il979h and 
star counts for larger radii from iKing et al.l (jl968t) . The 
star counts are characterized by a magnitude limit of B=19 
mag. We use the star catalog of AvdMlO to perform star 
counts in the central regions (i? < 20") with respect to the 
kinematic center described above. We use stars brighter 
than 19.5 mag and an adaptive kernel density estimator for 



the star counts. The magnitude cut is applied to limit the 
incompleteness. Our profile center is stable for this mag- 
nitude cut-off (fainter magnitude cut-offs cause the density 
center to shift towards the Anderson center). We adjust our 
star count profile to the Meylan (1987) and Trager (1995) 
profiles at larger radii, as taken from NGB08. In NGB08 a 
surface brightness profile was obtained with integrated light 
from HST-ACS data within the cent ral 40" with respect 
to the ir ce nter (details explained in iNovola k, Gebhardtl 
(I2006h and INovola fc Gebhardti (I2007D 1. We use their data 
from radii larger than 20", i.e. the data from 20" to 40" 
comes from integrated light measurements and the inner 
20" comes from star counts. 



2.3. Kinematical Data 

NIO obtain kinematics in the central region of w Centauri 
using integral field spectroscopy. They measure the veloc- 
ity dispersion from integrated light using VLT-FLAMES 
with a spectral resolution of ~'10,000 in the Ca-triplet 
wavelength. They tile around the two proposed centers by 
NGB08 and AvdMlO with eight pointings. We also use 
Gemini-GMOS data which NGB08 obtain with integrated 
light using the same approach as for the VLT-FLAMES 
data. We use the integrated light velocity dispersion with 
respect to the kinematic center as presented in Table 1 of 
NIO. 

vdV06 collect individual velocity measurements at 
larger r adii from four diffe r ent sources (iSuntzeff fc KraftI 
(I1996D : iMavor et al.l (|1997D : iReiins et al.l (|2006[ ): Xie and 
Gebhardt (private communication)). Almost all of the 
above authors measure the velocities of luminous (giants) 
stars. vdV06 perform many tests such as cluster member- 
ship, excluding velocities with large uncertainties and also 
corrections for perspective rotation, in order to pick only 
suitable velocities for dynamical modeling. They bin the 
measurements and obtain the velocity moments in a set of 
apertures in the plane of the sky. We use the velocity dis- 
persions presented in their Table 4 for comparison with our 
simulations. 

Proper motions are very useful to better constrain the 
internal dynamics of star clusters, in particular the de- 
gree of anisotropy. In addition to ground-based data from 
vdV06, HST proper motions are available from AvdMlO. 
These authors use isolated stars brighter than the apparent 
magnitude 21 in their high quality sample for proper mo- 
tions. In total, they have about 72,000 stars in two fields: 
one on the cluster center and one positioned adjacent to 
the first field along the major axis. The central field covers 
the central 147" in radius, and the major axis field covers 
radii between about 100" to 347". They use 25,167 stars at 
R < 71". 7, aiming at having the complete position angle 
coverage in order to calculate average kinematical quan- 
tities over circular annuli. However, they stress that the 
whole data set is usable but is excluded in their main study 
because of sparse position angle coverage at larger radii. 

To compare with the N-body models in this work, we 
use the proper motions on the minor and major axes avail- 
able in Table 4 of AvdMlO transformed with respect to 
the kinematic center in NIO. We measure the proper mo- 
tion dispersion along each axis using a maximum likelihood 
technique in radial bin s taking uncertainties into account 
(jPrvor fc Mevlanlll993D . 
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Table 1. biitial parameters, main assumptions and observed properties used in our models. 



Property 


Symbol 


Values 


Num. of stars 




5 X 10* 


Structure Model 




King (1966) 


Initial concentration 


c 


0.3 - 0.8 


Initial half-light radius 


initial r^p 


11.6 - 14.0 pc 


Initial mass function 


Kroupa 2001 


0.1 - 100 Mq 


Tidal field 




none 


Primordial binaries 




none 


Primordial mass segregation 




none 


Mean metallicity 


[Fe/H\ 


-1.62 (Harris 1996) 



" see Section 3 for scaling description. 



Throughout this work, we assume a heliocentric dis- 
tance of 4.8 ± 0.3 kpc for lo Centauri (vdV06). Therefore, 
1 pc corresponds to 42.97" in our simulations. 

3. A^-body Modeling Method 

We started running simulations on GRAPE spe cial pur- 
pose computers at ESO using the NB0DY4 code (|A arse t3 
|1999( ) in order to model the star cluster uj Centauri. It be- 
came possible in the middle of this project, however, to 
take advantage of the recently installed GPU cluster of the 
University of Queensland which speeds up the simulations 
by a factor of about 10, allowing us to probe a larger ini- 
tial parameter space considerably cheaper in computational 
time. All the results of our models for w Centauri are ba sed 
on simulations using the NB0DY6 code H arse tbl [200l on 
the GPU cluster. 

We fo llow t he method described in iBaumgardt et all 
(|2003al lbl. 120051 ) to model the dynamical evolution of w 
Centauri. We set up our model clusters following a spheri- 
cal isotropic King model (|Kinelll966D in vir ial equilibrium . 
The initial stellar masses are drawn from a iKroupal (|200l[ ) 
initial mass function (IMF) with lower and upper mass lim- 
its of 0.1 and 100 Mq, respectively. Such a choice of IMF 
is supported in Section 4 since it reproduces well the M/L 
profile (see Section 4.2), consistent with other independent 
studies (vdV06). Primordial binaries are not included in 
our simulations. This is justified since lo Centauri has a 
very high velocity dispersion, which reduces the contribu- 
tion of a reasonable small fraction of primordial binaries. 
In addition, including primordial binaries is computation- 
ally very time consuming and we will investigate it in a 
separate project. We also do not include initial mass seg- 
regation. Furthermore, we do not consider the tidal field 
of the Galaxy since we are interested in the very central 
part of the cluster. Neglecting the infiuence of the tidal 
field, the evolution of our modeled clusters is driven mainly 
by two-body relaxation and stellar evolution. We note that 
mass loss has little effect on the current velocity dispersion 
profile if the mass loss occurred early in the evolution. In 
our simu lations, stellar evolution is modeled according to 
iHurlev et al. (2003). We assu me [Fe /H] = -1.62 dex as 
the mean cluster metallicity ()Harrisl Il996). The assumed 
neutron star and black hole retention fraction is set to 10% 
for both no-IMBH models and models with an IMBH. The 
initial parameters in our simulations are summarized in 
Table [fl 

Centauri has a mass of ~ 2.5 x lO^M© and therefore 
~ 5 X 10^ stars. Since direct N-body simulations can at 



the moment handle only cluste rs with up to 10^ stars (see 
e.g.^Hasani Zonoozi et al.ll2010t) . we perform simulations of 
smaller-N clusters (but more extended in size) and scale 
our results up to uj Centauri such as to have the same re- 
laxation time and size as our observed cluster after 12 Gyr 
of evolution. The relaxation time of a cluster with mass M 
and half-mass radius is given bv lSpitzed (ITqSTI) as: 



Tr^ = 0.138- 



'M r 



3/2 



m) VG ln(7A^) 



(1) 



where (m) is the mean mass of the stars in the cluster, N is 
the number of stars, and 7 is a factor in the Coulomb loga- 
rithm, approxima.tely e qual to 0.02 for multi-mass clusters 
()Giersz fc Heggiell 19961 ) . The scaling factor for the positions 
is given by 



'^hpocen 



(2) 



where rhpocen is the projected half -light radius of uj 
Centauri (5.83 pc at 4.8 kpc distance, iHarrisI [79960 where 
the integrated cluster light reaches half its maximum value 
and rfipsim is the projected half-light radius of model clus- 
ters. In order to have the same relaxation time as the ob- 
served cluster, we have to scale up the mass of our clusters 
to a mass Mocen which satisfies the following equation: 



Ms^rn /ln(7Mocen)A^^'^ _ '''hp 



Mo 



'^hpsim 



(3) 



Here Mocen is the bound mass of the models at the end of 
our simulations (T = 12 Gyr). Since the size and mass of 
the model clusters are re-scaled, we also scale the velocities 
of the stars by a factor: 



scale 



\^'^ /M 



T^hpocen 



M,. 



(4) 



where subscripts "ocen" and "sim" denote the actual values 
for U J Centauri and those in the simulations, respectively 
(see IBaumgardt et al.ll2003bl ). In eq. (4), the first factor 
is needed due to the reduction of distances between stars 
while the second factor takes care of the increase in cluster 
mass when scaling our models to w Centauri. Due to eqs. (3) 
and (4) models starting with larger initial rhpsim will end 
with a higher cluster mass and therefore higher velocities 
(see Section 4). 

We perform all simulations with 5 x 10'' stars, which 
before scaling gives us an initial cluster mass of ~ 3.5 x 
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1O^M0 and a final cluster mass of ^ 2.0 x IQ'^Mq. Usually, 
the bound stars are about 99% of the initial number of 
stars. 

In order to find the initial conditions of w Centauri that 
lead to the present-day observed kinematics and surface 
brightness profile, we set up clusters with different initial 
half-mass radii, different concentration parameters, defined 
as log(^) where rt is the tidal radius and Tc is the core 
radius, and different IMBH masses and evolve them with 
NB0DY6 up to an age of 12 Gyr. Then, we estimate how 
closely each model cluster reproduce the observed profiles 
using a test. We calculate a value for the surface 
brightness, the radial velocity, and the proper motion dis- 
persion profiles of each model to compare with the data 
using 



X 



N 

E 



(5) 



M.i, Vi, AAii and AVi are the model and data points and 
their relevant uncertainties. N is the number of data points 
for the radial velocity dispersion, the proper motion disper- 
sion and the surface brightness profiles. We calculate the 
absolute values for all the maps in the next section. 
We use values only to quantify judgments on different 
model profiles in comparison with observed ones. We cal- 
culate the values for all models within the inner 400" 



10 pc) since we do not consider tidal fields and this in 
turn affects the number of stars at larger radii. We aim to 
simultaneously reproduce the observed velocity dispersions 
and surface brightness profiles. Therefore, we apply the fol- 
lowing relation to obtain the reduced combined values 
of radial velocity, and proper motion dispersion and surface 
brightness for each model: 



X 



X^ 

A,rv 



Apm 



xlb 



Nr 



(6) 



"rv", "pm" and "sb" stand for radial velocity dispersion, 
proper motion dispersion and surface brightness, respec- 
tively. Nrv, Np^, and Nsb are 38, 218, and 32 over 400". 
We note that the large difference between the number of 
"rv" , "pm" and "sb" data points will cause a non-smooth 
total x^ space between different models and consequently 
results in small absolute differences, as we see in Section 
4.2. 

We vary the initial parameters to compare the result- 
ing profiles with observations of w Centauri. We vary the 
IMBH mass between 0% to 3% of the model stellar mass. In 
total we compute more than 100 models to find the initial 
conditions which reproduce the observations of ut Centauri 
best. The final quantities (e.g. number density and velocity) 
are calculated by adding up three model clusters with the 
same physical initial conditions but different random num- 
ber seeds and also superimposing 10 snapshots in the case 
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Fig. 1. Radial velocity dispersion profiles as a function of 
radius for models with an IMBH mass of 2% of the stellar 
mass with the same initial concentration factor. Red points 
are the observed data taken from different sources (section 
2.3). Black, magenta, light blue and green are models with 
initial half-light radius (rhp) of 11.6, 12.2, 12.8 and 13.4 
pc, respectively. A typical uncertainty for models at 5" is 
indicated in blue. Higher initial Vhp models end up with a 
higher velocity dispersion across all radii but with roughly 
the same overall shape because of their higher total cluster 
mass as described in section 4. 
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Fig. 2. Radial velocity dispersion profiles as a function of 
radius for models with an IMBH mass 2% of the stellar mass 
with the same initial r^p- Observational data are shown as 
red points, and a sample uncertainty is shown in blue as 
in Figure [TJ Black, magenta, light blue, green and orange 
points are models with initial c — 0.4 to 0.8. All models 
have the same dispersion outside the core radius (~ 100") 
while clusters with higher concentration show higher veloc- 
ity dispersion towards the center, as explained in section 
4. 



6 



B. Jalali et al.: A Dynamical N-body Model for the Central Region of uj Centauri 



of no-IMBH models and 20 snapshots in the case of IMBH 
models to improve the statistics. The snapshots start at 11 
Gyr with a 50 Myr step. 

3.1. Model Density Profile 

For each model we calculate the surface brightness and 
kinematic profiles including radial velocity and proper mo- 
tions, using similar magnitude limits as observers use for 
u Centauri. In order to determine the surface brightness 
profile, the density center o f our model clusters is d eter- 
mined using the method of ICasertano fc HutI (|1985[ ). We 
then count the number of bound stars in two-dimensional 
circular annuli around the density center. We use the infi- 
nite projection method of iMashchenko fc Silld (|2005f ) and 
average each quantity over all (infinite) orientations for each 
bin based on geometric arguments. Infinite projection gives 
significantly better statistics over using only a finite number 
of projections, especially in the inner cluster parts. We con- 
vert the star counts per parsec squared to numbers per arc- 
second squared using the assumed distance to w Centauri. 
We bin the stars around the (density) center in 20 annuli 
of equal logarithmic width between 2.0" and 1000". 

We separately calculate the surface number density of 
all stars (including dark remnants) and only bright stars. 
In order to compare the models with the observations, we 
only consider stars brighter than \—22 magnitude. We con- 
vert model bolometric luminosities to V-ba nd luminosities 
assum ing the stellar temperature model of lEggleton et al] 
()1989f ). Using the distance modulus of u Centauri, we ob- 
tain V-band magnitudes that can be directly compared to 
the measurements (e.g. magnitude cuts). AvdMlO measure 
the observed surface brightness profile using HST multi- 
epoch data with the magnitude cut of B~22 mag (section 
3.2 and their Figure 1). Our magnitude cutoff of V= 22 mag- 
nitude is al so consistent with th e combined data in MevlanI 
(|1987^ and lTraeer et all (|1995D . The adopted cutoff in our 
simulations is also applied to maximize the number of stars 
in each bin for better number statistics. We then convert 
the star counts to magnitude per arcsecond squared. We 
match our surface brightness profile with the observed one 
by shifting it by an additional zero point. We calculate the 
zero point for each model profile with a minimization in 
order to compare our profile directly to the observed surface 
brightness profile. 

3.2. Model Kinematic Profile 

We calculate the velocity dispersion profile of our model 
clusters again using all stars and using only bright stars 
similar to how the surface brightness profile was calculated. 
First, we determine the velocity dispersion using all stars 
including compact remnants. Second, we use only stars 
brighter than a certain magnitude limit, adopted to be the 
same as the observational one. In the case of radial ve- 
locities, the observed kinematical data within the central 
30" are obtained using integrated light (IFU data), so we 
similarly measure in the models the luminosity weighted 
velocity dispersions. We avoid very bright stars in the cen- 
tral region as in the observations, in order to minimize shot 
noise effects. Therefore, we only consider stars fainter than 
V-magnitude of 15 (similar to the observational cut off) 
inside 30". For radial velocities at larger radii, we consider 



stars brighter than an apparent magnitude of V=18 magni- 
tude since the observed data were mainly obtained from in- 
dividual giant stars. For proper motions, we consider stars 
within the magnitude range of 18 < V < 22 magnitude 
which is found suitable for proper motion measurements 
in AvdMlO. We assign magnitude weights to the velocities 
of stars within this magnitude range following Table 4 of 
AvdMlO to measure the proper motion velocities. 

4. Results 

In order to find a model which has simultaneously a good 
fit of the observed surface brightness and velocity disper- 
sion profiles of UJ Centauri, we run a grid of models with 
different initial conditions. We describe below some general 
phenomena in order to illustrate the effect of each initial 
parameter on the profiles of the evolved clusters. We use 
as an example the IMBH model with 2% mass of the to- 
tal cluster mass. We first consider variations of the initial 
rhp while the cluster concentration (c) is fixed. The model 
cluster final mass is a free parameter as discussed for eq. 
(3). Since the final Vhp is fixed t o the oj Centa uri one (4.18 
arcmin in the 2003 version of [Harris! (fTgOeh l. increasing 
the initial radius will produce a more massive cluster af- 
ter scaling, since we scale such as to keep the relaxation 
time constant (see eq. 3). Hence, at a fixed initial cluster 
concentration, by increasing the initial r^p the whole radial 
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Fig. 3. Upper panel: V-band surface brightness profile as a 
function of radius. The red points are the observed V-band 
surface brightness relative to the kinematic center in NIO. 
The black points are the best no-IMBH model. This model 
fits the surface brightness data relatively well inside 400", 
except inside the central 10", where it ends up below the 
observations. The model points are connected inside 400" 
over which we calculate the values. Lower panel: residual 
of our model with respect to the observed profile. 
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Fig. 4. Upper panel: Radial velocity dispersion profile vs. 
radius. The red points are the observed velocity dispersion 
relative to the kinematic center, taken from NIO. The ve- 
locity dispersion of the best-fit no-IMBH model is shown in 
black. The model obviously does not fit the data around the 
center (for details see section 4.1). Lower panel: the relative 
difference of our model and the observed profile. 
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Fig. 5. The proper motion velocity dispersion profile vs. ra- 
dius. The best-fit no-IMBH model is shown in black. Shaded 
magenta and green are the observed proper motions for ma- 
jor and minor axes taken from vdMAlO but with respect to 
the kinematic center. The Lower panel shows the residuals 
and displays the major axis data only, for clarity. 



velocity dispersion profile will scale up, as can be seen in 
Fig HI In this figure, at a fixed concentration c — 0.4 one 
can see that a cluster which starts with an initial r/jp of 
11.6 pc has a much lower velocity dispersion profile than 
the observed one at almost all radii. In contrast, the cluster 
with the same concentration but higher initial r^^ of 13.4 
pc has a higher dispersion profile than the observed one, 
while the general shape of the profiles usually follows the 
same pattern. This also shows that relaxation is not very 
important for a; Centauri. 

We now discuss the effect of varying the initial concen- 
tration (c) on the final kinematic profiles after 12 Gyr of 
evolution. Clusters with higher concentrations have more 
mass in the central regions and therefore a higher cen- 
tral velocity dispersion. In Figure [2j we show this effect 
by presenting one family of models with fixed initial Vh-p 
(12.2 pc) but varying concentrations (see the color code in 
the caption). In this example, all clusters in the family of 
Thp = 12.2 pc have almost the same velocity dispersion at 
large radii, but different central velocity dispersions as a 
function of their concentrations. 

4.1. No-IMBH Models 

We first run a sparse grid of models between 0.5 < c < 1.5 
and 10.0 < r^p < 15.0 pc in order to identify the best 
fitting model. We find that models with initial c ^ 0.8 and 
rfip ^ 12 pc give the best fit. We produce a finer grid of 
models between 0.6 < c < 1.0 and 12.2 < r^p < 14.0 



pc for the no black hole case. As explained in section 3, 
we calculate the values of velocity dispersion, proper 
motion dispersion and surface brightness profiles to choose 
the best-fitting no-IMBH model. We find that the model 
with initial c = 0.8 and initial r^p = 12.8 pc produces the 
best fit to the data of uj Centauri. 

Figure [3] shows the surface brightness profile of bright 
stars with V < 22 magnitude for the model with c = 0.8 
and Thp = 12.8 pc and no central black hole. In this model, 
the surface brightness of bright stars reproduces the obser- 
vational data well but predicts slightly lower values in the 
central 10" where NGB08 see a shallow cusp, which can be 
interpreted as an evidence of an IMBH. Fig. 2] shows the 
radial velocity dispersion profile for the above no-IMBH 
model. Except for the innermost part, the model agrees 
relatively well with the data within the uncertainties but 
in the central 100" the no-IMBH model lies significantly be- 
low the observational data. The proper motion dispersion 
profile for the no-IMBH best-fit model is shown in Fig. [S] 

We note that whether the no-IMBH model does a 
good or a poor fit depends on the adopted center since 
in vdMAlO there is no shallow cusp in the observed sur- 
face brightness profile. Further, the observed proper motion 
data around the proposed center by AvdMlO would be bet- 
ter fitted by the no-IMBH models although the central 10" 
values are still slightly above the best no-IMBH models. 

Models starting with higher r^p than 12.8 pc cannot 
match the radial (and proper motion) velocity dispersion 
profile beyond ~ 50", they always lie higher than the 
data at larger radii when we attempted to match the in- 
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Fig. 6. Velocity dispersion profile with different central 
black hole masses. The red points are the observed dis- 
persion. The black line is the best no-IMBH model which 
obviously cannot reproduce the data in the inner parts. The 
blue, green and yellow lines represent models with black 
hole masses of 1%, 2% and 3% of the cluster total mass. A 
typical uncertainty at 5" is shown in blue on top of models 
for visual clarity. The green line is the velocity dispersion 
profile for the best-fit black hole model, containing 2% of 
the initial cluster stellar mass (section 4.2). 



ner region. In other words, higher initial r^p, meaning 
higher model cluster mass, increase the whole model profile. 
Therefore, models with an initial r^p higher than 12.8 pc are 
ruled out. Similarly, models with 12.2 pc initial rhp do not 
reproduce the observations either. However, the majority 
of the data points used to compute the radial velocity dis- 
persion is at radii larger than ~ 100", while for the proper 
motion dispersion it is at radii smaller than ^ 100". This 
can lead to a slightly different best-fit model depending on 
the radial range over which the model is evaluated. Given 
eq. (3), we use the data at large radii to choose the best 
cluster mass scaHng. For instance, the model with c = 0.8 
and rhp = 12.8 pc better fits the proper motions than the 
c = 0.8 and r^p = 12.2 pc which better fits the radial veloc- 
ity dispersion. The latter provides a good fit to the data, 
but the total is slightly larger than for the former since 
the number of proper motion data points is much larger 
(218 points) than the number of data points for radial ve- 
locity data (38 points), giving a higher weight to the proper 
motion fit. 

4.2. IMBH Models 

Since models without an IMBH could not represent the 
data well, we run models including central black holes of 
various masses with the hope of improving the fit to the 
data. In this set of simulations, we start from isotropic 



Fig. 7. Surface brightness profile of models with different 
central black hole masses. Symbols and colors are the same 
as Fig m The red points are the observed profile (section 
2). The black line is the best no-IMBH model which falls 
slightly below the observed data. The blue, green and yellow 
lines represent models with black hole masses of 1%, 2% and 
3% of the cluster total mass. 

King model conditions with concentrations in the range 
0.3 < c < 0.8 and initial rhp in the range of 12.2 < Vhp < 
14.0 pc. Apart from the central black hole mass, all the 
other parameters such as NS retention fraction and stel- 
lar mass range are the same as for the no-IMBH models. 
Furthermore, since we find that the IMBH model that re- 
produces the observed profiles best has lower concentra- 
tion and higher rhp than the no-IMBH models, we explore 
a larger parameter space in order to determine the best 
set initial conditions leading to the present-day observa- 
tions. For each initial configuration, we run the simulations 
as described in Section 3. For models with a central black 
hole, we calculate the kinematic and surface brightness pro- 
files as described in the Sections 3.1 and 3.2. Hence, all the 
model-data comparisons are done following the same mag- 
nitude and radial cut-offs as for the no-IMBH models. In 
the case of IMBH models, the IMBH moves in the core of 
the cluster during the simulations, and we center the final 
model clusters on the central black hole. In our simula- 
tions, IMBH models contain black hole masses of 1%, 2% 
and 3% of the final stellar mass of the clust er. For reference, 
a 1% IMBH would lie slightly above the iMaeorrian et al.l 
(11,998) relation. If we adopt a mass of ~ 2.5 x lO^Mg for lu 
Centauri (vdV06), the IMBH masses we use are equivalent 
to black hole masses of about 2.5x lO^M©, 5.0x lO'^M© and 
7.5 X 1O^M0 respectively. However, the exact masses of our 
models depend on the initial rhp due to the scaling of radii 
(see eq. 3), and will slightly differ from the above values. We 
calculate and analyze the maps for models including the 
above IMBH masses. We find the best-fit models for each 
grid including a 1%, 2% and 3% IMBH, applying the same 
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Fig. 8. From left to right, the map for the radial velocity, the proper motion dispersion and the surface brightness 
profiles, for models with an IMBH mass of 2% of the model cluster mass, over the grid of initial parameters space. The 
numbers in black are the absolute over 400" radius. The absolute values inside 40" are shown in red to visualize 
the goodness of the fit in the central region for each model. 



methodology that we use to find the best no-IMBH model. 
The radial velocity dispersion and surface brightness pro- 
files of the best-fit models of different IMBH masses and 
also the ones for the best-fit no-IMBH model are shown 
in Figure IH and Figure [T] The best-fit IMBH model is the 
one containing 2% of the stellar cluster mass starting ini- 
tially with c = 0.5 and rtp = 13.4 pc. As can be seen from 
Figureini the best-fit IMBH models with 1% and 3% of the 
cluster stellar mass do not fit the observed radial velocity 
dispersion well since they have lower or higher values than 
the observed velocity dispersion profile in the central part. 
A summary of best-fit models among all models including 
no-IMBH and IMBH are tabulated in Table H The chosen 
best-fit 2% mass IMBH model is highlighted with boldface. 

The left panel in Figure |8] shows the radial velocity dis- 
persion map for models containing an IMBH mass of 
2% of the cluster stellar mass. Models starting with an ini- 
tial rhp of 12.2 (14.0) pc always lie lower (higher) than the 
data at large radii because of the less (more) massive final 
cluster mass (see eq. 3). In contrast, models starting with 
Thp of 12.8 and 13.4 fit the data much better at larger radii 
due to better cluster mass scaling when relaxation time is 
the same as for uj Centauri. 

In addition to radial velocities, we calculate the proper 
motion dispersion (perpendicular component to the radial 
velocity in our models) with the weights and magnitude cut- 
offs explained in section 2.3 and 3.2. The map for the 
proper motion dispersion for the 2% IMBH models is shown 
in the middle panel of Fig. [S] The right panel in Figure [S] 
depicts the x^ values of surface brightness for the mod- 
els with an IMBH mass of 2% of the stellar mass. Models 
starting with high initial concentrations such as c = 0.7 and 
0.8 generally show a steeper cusp in the central part than 
the observation and do not describe the observed surface 
brightness well, they have higher x^ values. 

A larger core gives a better fit at larger radii (around 
80"), though these models have a very poor fit in the cen- 
tral part due to a very steep rise. In addition, it should 
be noted that the number of data points at large radii is 
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Fig. 9. The final reduced, combined x^ map for models 
with an IMBH with a mass of 2% of the model cluster 
mass. This map is calculated based on equation 6, using 
the radial velocity, the proper motion dispersion and the 
surface brightness x^ maps in Fig [8] The best- fit model is 
marked in filled green. 

higher than at the center. Furthermore, the observational 
uncertainties are smaller at larger radii than in the central 
region, therefore, the x^ value can be small if a model has 
a good fit at large radii but an unsatisfactory fit in the cen- 
tral part. Thus, we choose models with smaller (c) such as 
0.4 and 0.5 as our best-fit models for the surface brightness 
data. The fact that smaller core models such as c = 0.4 
and 0.5 do not provide a better fit at larger radii in surface 
brightness is related to the initial King profile, which alone 
might not be sufficient to fit the density profile and is in- 
dependent of having a central IMBH. In a future paper we 
intend to modify the initial density profile so that we take 
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30 



Fig. 10. Upper panel: V-band surface brightness profile as 
a function of radius for the best-fit model with an IMBH 
mass of 2% of stellar mass. The symbols are as in Fig. [TT] 
Lower panel: residual of our model and the observed profile. 



this effect into account, for instance by combining a King 
profile with a Sersic one. 

Considering the above results, the model with initial 
c = 0.5 and rtp = 13.4 pc fits better than the other models. 
However, it still has a slightly steeper surface brightness in 
the central region than the observed one. We show the final 
reduced combined values (see eq. 6) for models with 
a 2% IMBH mass in Figure The IMBH model starting 
with initial c — 0.5 and rhp — 13.4 pc provides the best- 
fit to the observations. The profiles for the best-fit IMBH 
model are shown in Figures [TUl [11] and [T^l As can be 
seen in Figure [T21 a model with an IMBH mass of 2% of 
the final cluster stellar mass shows higher central velocity 
dispersions than the observed proper motions. However, as 
mentioned before, the best-fit no- IMBH model might fit the 
proper motions better with respect to the AvdMlO center. 

Furthermore, it is interesting to note that the radial ver- 
sus transversal components of the velocities for our evolved 
clusters are the same. Thus, an initially isotropic cluster 
stays isotropic throughout its evolution even if a central 
IMBH is present. 

In Fig. [131 we show the stellar V-band mass-to-light ra- 
tio profiles in our simulations for the best-fit models with 
an IMBH, and the model without an IMBH. Since the ab- 
solute value of M/L depends on the age of the stars, we 
calculated it using the same 20 snapshots, i.e. from 11 to 
12 Gyr, for ah the models in Fig. [H The M/L ratio is the 
same for all models at large radii since the initial mass func- 
tion is the same, and, the effect of an IMBH is negligible 
in the outer region. The M/L ratio seems to be constant in 
the central region and is increasing outward starting from 
the core radius (~ 100"). It is about 10% larger than its 
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Fig. 11. Upper panel: Velocity dispersion profile vs. radius 
in arcsecond. The red points are the observed data points 
relative to the kinematic center, taken from NIO. The veloc- 
ities for the best-fit model containing an IMBH mass of 2% 
of stellar mass are shown in black. Lower panel: the relative 
difference between our model and the observed profile. 



value at the core radius for all models due to the increase 
in number of low mass stars at larger radii. 

As mentioned in Section 3, setting the relaxation time 
and fixing the final present-day rhp to the one of a; Centauri, 
we obtain the final (present-day) total mass of the model 
cluster. We derive a cluster mass of (2.6 ± 0.1) x lO^M© if 
we use the best-fit model containing an IMBH of 2% of the 
cluster mass. This is in very good agreement with the mass 
of cj Centauri as determined by vdV06, (2.5±0.3) x 10*5 Mg. 



5. Discussion and Conclusions 

We have created a large set of evolutionary N-body mod- 
els for Lo Centauri assuming that the cluster started with 
a Kroupa IMF, with the aim of reproducing its observed 
properties. The main goal is to check models using standard 
assumptions as a first attempt, and to see to what extent 
we reproduce the observations. In particular, we examine 
whether we can explain the newly acquired observations for 
the central velocity dispersion profile with the presence of 
a central IMBH or whether a model without an IMBH is 
also consistent with the observations. 

Following the method applied in Baumgardt et al. 
(2003a, b), we calculate models starting with spherical 
isotropic King conditions (King 1966) with different ini- 
tial parameters. We do not include the tidal field of the 
Galaxy and primordial binaries in our models. Since we can- 
not simulate a star cluster of the size of a; Centauri by direct 
N-body simulations, we start with more extended clusters 
containing fewer number of stars than io Centauri and scale 
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Table 2. Initial parameters and the values for best-fit models containing different IMBH mass fractions. 
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Model 


c 


rhp (pc) 


absolute xtv 


absolute Xpm 


absolute xtb 


reduced Xtotail 


reduced xLtai 


40'X 


no-IMBH 


0.8 


12.8 


135.61 


294.80 


537.31 


3.36 


3.06 




1% IMBH 


0.7 


12.8 


99.91 


345.52 


313.21 


2.63 


3.80 




2% IMBH 


0.5 


13.4 


62.41 


271.51 


277.95 


3.09 


2.58 




3% IMBH 


0.4 


13.4 


58.93 


325.88 


654.72 


3.61 


5.01 





" This column is obtained by dividing the sum of absolute on the total number of data points for all three observed profiles 
' The same as the previous column but for the inner 40" region 
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Fig. 12. The proper motion dispersion profile for our best- 
fit IMBH model is shown in black. Shaded magenta and 
green are the observed proper motions for major and minor 
axes taken from vdMAlO but with respect to the kinematic 
center. The model proper motion dispersions for a cluster 
containing an IMBH mass of 2% of the stellar mass has 
a higher velocity dispersions than the observed one in the 
central 10". The residual is illustrated only for major axis 
data for clarity. 



our model clusters to the observed cluster such that the re- 
laxation time is constant. As described in detail in Section 
3, we measure physical quantities from our models such 
as velocity dispersion and surface brightness following the 
methods in observational studies, with the same magnitude 
cut-offs and luminosity weights. Using such careful magni- 
tude weights and radial cut-offs in our models make the 
data-model comparisons and consequently the drawn con- 
clusions more reliable. We use values to compare the 
profiles of model clusters after 12 Gyr of evolution with the 
observed ones. In Section 4, we present a grid of models 
for clusters containing an IMBH mass of 2% of the cluster 
stellar mass. 

We show that the best-fit IMBH model, containing a 
5 X 10'* M© black hole, matches the data presented in Noyola 
et al. (2010) very well. In particular, we reproduce the ob- 
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Fig. 13. V-band stellar M/L ratio as a function of distance 
to cluster center for the best-fit models normalized by M/L 
at half-light radius. The M/L ratio increases slightly out- 
ward from the core radius due to an increase in the fraction 
of low mass stars. 



served rise in the central velocity dispersion as an indicator 
for the presence of an IMBH. We stress that relying on 
the profiles relative to the kinematic center in Noyola et 
al. (2010) makes it impossible to consistently fit the radial 
velocity dispersion over all radii without a central IMBH. 
Furthermore, we show that M/L is increasing from the core 
radius towards large radii for all models independent of the 
presence of a central black hole. We do not claim that our 
best-fit model is unique. However, we examine more than 
100 models to be confident about the initial parameters and 
the final chosen best-fit model. 

A number of further details can be investigated or im- 
proved as a next step. Among the interesting issues are the 
study of rotation at different radii, especially in the cen- 
tral region. In principle, such a study can be done using 
an axisymmetric model which is close r to uj Centauri, as 
an ex ception in Galactic star clusters. iFiestas fc SpurzemI 
(|2010D investigated the evolution of rotating dense stellar 
systems containing massive black holes. Exploring rota- 
tion effects could help to better understand the observed 
discrepancy between proper motion and radial velocity 
dispersion. In addition, as described in previous studies 
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(|King fc Anderson! 120021: iMevlanI IToSTt ) . oj Centauri might 
not be well fitted by a King profile alone. We intend to 
investigate this in more detail in a future paper studying 
different initial configurations such as Sersic and double 
King profiles. Furthermore, to improve the model surface 
brightness profile at large radii, around the tidal radius, the 
tidal field of the Galaxy should be taken into account. In 
order to more tightly constrain the initial parameters for w 
Centauri the above studies as well as trying different IMF 
distributions are necessary in addition to investigate a finer 
grid, though they will be computationally expensive. 
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